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Abstract 

We derive expressions for the lossy boundary-scattering contribution to the linewidth of surface 
electronic states confined with atomic corrals and island resonators. Correcting experimentally 
measured linewidths for these contributions along with thermal and intrumental broadening en- 
ables intrinsic many-body lifetimes due to electron-electron and electron-phonon scattering to be 
determined. In small resonators lossy-scattering dominates linewidths whilst different scaling of 
widths and separations cause levels to merge in large resonators. Our results enable the design of 
resonators suitable for lifetime studies. 

PACS numbers: 73.20.At,79.60.Jv,72.10.Fk,73.21.-b 



I. INTRODUCTION 



The dynamics of hot electrons and holes at surfaces have attracted considerable attention 
in recent years, not least due to their importance in photochemistry and charge-transfer pro- 
cesses such as electronically-induced adsorbate reactions. A particular focus has been the 
study of quasiparticle lifetimes associated with excitations in the band of surface-localised 
electronic states found at the (111) surfaces of the noble metals, which have become a test- 
ing ground for new theoretical, computational and experimental methods aimed at devel- 
oping a deeper fundamental understanding of quasiparticle dynamics. 1 Theoretically, much 
progress has been made towards a quantitative account of inelastic electron-electron^ (e-e) 
and electron-phonon^ (e-p) interactions of these states, whilst techniques such as scanning 
tunneling microscopy^ (STM) and spectroscopy 6 - (STS), and photoelectron spectroscopy^ 
have been used to construct a steadily-increasing database of experimentally-determined 
lifetimes. The lifetimes of excitations with energies spanning just a few tenths of an eV have 
been shown to reflect a wealth of key surface physics; intra and interband scattering pro- 
cesses; spatially-dependent and <i-electron screening; defect-scattering processes; electron- 
phonon interactions with Rayleigh (surface) and bulk phonons; and their temperature de- 
pendencies. 

Lifetime studies using the STM can be divided into two approaches: those which exploit 
the phase coherence of quantum interference patterns^ and those based upon lineshape 
analysis 6 . The lineshape method uses spectroscopic measurements of the differential con- 
ductivity dl / dV at a fixed position above the surface, and relies upon the presence of spectral 
structure that contains a signature of the quasiparticle lifetime. On pristine surfaces the 
only such feature is the band edge onset, so that only the lifetime of excitations at this en- 
ergy are accessible. However, electron confinement to natural or artificial nanoscale electron 
resonators, such as those shown in Figure ^ results in energy quantisation, inducing spectral 
structure in the form of a series of resonant levels at energies that can be controlled through 
hanges in the dimensions and geometry of the resonator. The quasiparticle lifetime is then 
reflected in the level widths*^, but additional contributions arise due to lossy boundary 
scattering that must be accounted for if lineshape analysis is to be used to determine the 
intrinsic quasiparticle lifetime. In this paper we provide analytic and numerical results for 
the linewidths of electrons confined by different nanoscale resonators, and discuss their use 
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in lifetime studies using STS. 

The starting point of our analyses is the equation satisfied by the single particle Green 
function G(r,r') n 

f--^+V-£+iIm£ I J G(r,r') = -5(r-r') (1) 

where r is a two-dimensional position vector, m* the effective mass and E the electron energy 
with respect to the surface state band minimum. Confinement is due to the potential V, 
which in general will remain unspecified. Only the scattering properties of V are relevant. 
We also include an inelastic potential or self energy ilmHi(E) to account for the effects 
of e-e and e-p scattering. This is related to the lifetime t\ associated with these processes 
by T\ = — ft/(2ImEi)ji The local density of surface states (LDOS) is obtained from the 
Green function as p(r; E) = — (2/7r)ImCr(r, r) (the 2 is for spin-degeneracy). We associate 
the LDOS with dl/dV in the usual way^. In section |H] we consider circular atomic corral 
resonators, and in section IIHI circular adatom and vacancy islands. Non-circular resonators 
are discussed in section IIVI Contributions to experimentally measured linewidths arising 
from thermal broadening and instrumental effects are described in section and in section 
I VII we illustrate how the results of our analyses may be used to identify appropriate resonator 
structures for lifetime studies. 




FIG. 1: Nanoscale electron resonators used in lifetime studies, (a) Pseudo three-dimensional 
representation of a constant current STM image of a rectangular corral approximately 9 nm x 
10 nm constructed from 28 Mn atoms on Ag(lll). In order to enhance the atomic structure the 
intensity is plotted on a logarithmic scale. The standing wave pattern inside is due to electron 
confinement^, (b) Greyscale representation of a constant-current STM image of a 240 nm x 240 nm 
region of a Ag(lll) surface showing numerous hexagonal monatomic vacancy islands. Surface state 
electrons are confined within these islands by strong scattering at the step edges of the vacancies. 
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II. ATOMIC CORRALS 



Atomic corrals are artificial structures in which individual adatoms are positioned with 
atomic scale precision into closed geometrie o 5 i 10 i 13 i 14 . Strong scattering of the surface state 
electrons by the surrounding adatoms leads to confinement. The "standard model" for 
describing these systems is the s-wave scattering modet^, in which the solution to (JTJ) is 
given b y 10 i 15 

G(r,r') = G (r,r')+J2 G o( r ^ R j) TJkG o(Rk,r') (2a) 

j,k 

T jh = t6 jk + t G (Rj, Ri)T lk . (2b) 

The sums in (j2J) are over the N identical adatoms at locations Rj, j = 0, 1, ..N—l, which are 
characterised by the scattering t-matrix t = (i/m*)(exp2i5 — 1), where 5 is the phaseshift. 
Go denotes the free-electron propagator— 

G,{ry) = -^H^{K\r-r'\) (3) 

where Hq is a Hankel function of the first kind^ 6 - and k = ^2m*(E — ilmTj{)/h 2 . In the 
case of a circular corral, and taking R., = S(cos$j, sin$j), dj = 2nj/N, when both r and r' 
are at the center of the corral 

and the LDOS at the corral center is 

where po = {m* / nh 2 ) (1 — (2 / n) arg/t) is the clean surface LDOS 6 - and 

S = -f^E^WRi-Rol) (6a) 
j'=i 

« ~ [nJ (kS)H^(kS) - C] . (6b) 



The last result follows using |Rj — Rq| = 25 sin •dj/ 2 and replacing the sum by an integral 
over -n/N < ft < 2n - tt/N. J is a Bessel function 16 and C = 1 + (2z*/7r)(ln(/«/4//) + 7 - 1), 
with 7 = 0.57721 56649... (Euler's constant). We have introduced ji = N/2ttS, the linear 



density of atoms making up the corral. Poles in the Green function which correspond to 
bound electron states that have amplitude at the center of the corral occur when t^ 1 = Q. 
Using the asymptotic forms^ for J , gives the condition as 



Q(k) = — 



C 



im*t s 1 



2h 2 



(7b) 



In the limit /x — > oo the solutions to ((Zj) are when k = n n = {nit — Tr/4)/S, corresponding to 
a series of resonances at energies 

h 2 (rm - 7r/4) 2 

En = 2m*S* (8) 
and with widths (full width at half maximum) 

T In = -2Im£/(£ n ) = /i/r ; (£ n ). (9) 

This is the "hard- wall" limit, where the surface state wavefunction vanishes at the radius of 
the corral.— The exact energies in this limit have the same form as (jSJ) but with (mr — 7r/4) 
(= 2.36, 5.50, 8.64 for n — 1, 2, 3 ) replaced by the n'th zero^ of the Bessel function J , j 0>n 
(= 2.40,5.52,8.65 for n = 1,2,3). This close agreement validates the use of asymptotic 
forms in obtaining (jZJ). 

In the hard-wall limit shows that the resonance widths are directly related to the 
inelastic e-e and e-p scattering rates at the energy of the level. However, in practice there 
are physical limits to how large the linear adatom density /i can become, with /i < where 
d nn is the nearest-neighbour separation of surface atoms. Consequently the confinement is 
always less effective than in the hard- wall limit, which results in increased level widths. 
Solving (J7J) to first order in this case yields level widths 

r„ = r In + r Cn (10) 

where 

r Cn ^^Reln[l-fi( Kn )]- (11) 
mo 

This is the contribution to the level widths due to lossy scattering by the corral adatoms. 
Figure El shows the typical behaviour of Tq, showing results for Fe corrals on Cu(lll) of 
various sizes and linear adatom densities. The accuracy of the approximations made in 




> 40 



£30 



CD 



60 



50 



10 








-0.4 -0.3 -0.2 -0.1 0.1 0.2 
E(eV) 



FIG. 2: Calculated level widths for Fe corrals on Cu(lll) using ImSi = 0, m* = 0.38 m e , im*t = 
1,Eq = —0.44 eV^ Squares and lozenges denote values found solving (J7J) by Newton-Raphson 
iteration for N = 48,5 = 71.3 A and N = 96, S = 142.6 A respectively. Circles denote values 
found by non-linear least-square fitting of a series of Lorenzians to the LDOS given by © for 
N = 48, S = 142.6 A. The lines are plots of (ti 2 k/ m* S)Reln[l - n(n)], k = ^/2m*(E - E )/H 2 for 
the respective corral parameters. 

obtaining the expression in (fTT|) are confirmed by the close agreement to be found between 
the level widths predicted by it and with those found by fitting the LDOS (0) with a series 
of Lorenzians or from iterative solution to (J7J). The general trend to be observed in the level 
widths is that they ncrease almost linearly with energy but are approximately halved by 
simultanously doubling the radius and the number of atoms in the boundary ring. 

III. ADATOM AND VACANCY ISLANDS 

In a second class of resonator confinement results from scattering at ascending or de- 
scending steps. These are adatom or vacancy islands . 17 i 18 i 19 i 20 In these systems the scatter- 
ing properties of the confining step are conveniently characterised by a reflection coefficient 



Assuming circular symmetry, with S now the island radius, the Green function can be 
expanded as 




(12) 



M 
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Substituting in to (JTJ) gives for g, g' < S 



2 d 2 G M dG M ( 2 2 n/ri\n 2m * xt >\ fio\ 

Q g 2 + + ( K e - M ) G M = -^-0 0(q-Q). (13) 

Solving by the direct method gives for g < g' < S 

tTTTTl* 

G m {q, q'\ E) = —— ttJm{kq)^! m {kq') (14) 

nr [A M — 1J 

where Jm is a Bessel function and ^> m{z) = H M \z) + AmHm\ z )i a linear combination of 
Hankel functions. The coefficient A M is chosen to ensure that ^ m satisfies the scattering 
boundary conditions at S. 

Identifying the energy levels from the poles in the Green function, we see from (JTljl that 
these arise whenever Am = 1- Only circularly symmetric states contribute to the local 
density of states at the center of the island ( Jm(0) = for M ^ 0) so that the energy levels 
that can be seen by STS measurements at the center of islands occur when Aq = 1. Using 
the asymptotic forms for the Hankel functions 

*o(«0) ~ .f^-le^- 71 ^ + A e- i(KS ~ n/A) ] 
V Ting L J 

K |y«(c-s) + Re -Mo-s)^ ; ( 15 ) 

that is, one can recognise within ty waves incident and reflected from the confining potential, 
enabling the coefficient Aq to be related to the planar reflection coefficient R of the step at 
the island edge. Using this relationship and equating A to 1 gives 

Re 2iKS e~ m/2 -1 = (16) 

as the condition for bound states visible in STS at the center of islands. Assuming confine- 
ment sufficient to give an identifiable series of resonant levels, and writing R=\R\ expi<ft R , 
predicts their energies as 

h 2 (nn + n/A-<p R /2) 2 
En = 2^ (17) 



and corresponding widths 



where 



... — Tin + Trh (18) 



h 2 2m*E n \n\R\ 
m* V h 2 S 



(19) 
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FIG. 3: Calculated level widths for circular Ag vacancy islands using ImS/ = 0, m* = 0.42, and 
Eq = —67 meV, for island radii 75 A(squares), 150 A(circles), and 300 A(lozenges). The lines are 



plots of — (ft?/m*S)^/2m*(E — Eq) /ft? In \R\ for the respective island radii. The inset shows the 
reflection coefficient used in the calculation (line), along with the experimental values from Ref. 

The "hard- wall" limit in this case is R — > —1, which ensures that the wavefunction (|15|) 
vanishes at S. In this limit (|17|) coincides with (jHJ) and r^ n = so that once again the 
resonance widths are directly related to the inelastic e-e and e-p scattering rates. In reality, 
scattering at real steps has been found to be lossy, so that \R\ < 1, and then (flUj) enables to be 
estimated the resulting additional contribution to the level width Tr. Figure El shows widths 
found by solving (fho|) for variously sized Ag vacancy islands using an energy-dependent step 
reflection coefficient fitted to experimentally determined values^, along with the widths 
expected using the approximate expression ()19jl. The agreement is again very good. The 
widths show similar behaviour to those in the atomic corrals, increasing approximately 
linearly with energy and inversely proportional to the island radius. 



IV. NON-CIRCULAR RESONATORS 

Equations (JTTj) and (flU|) enable the ready determination of the lossy scattering contribu- 
tion to the spectral linewidths of electron states in circular resonators. In terms of lifetime 
studies, circular resonators have two distinct advantages. Firstly, the high symmetry re- 
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FIG. 4: Calculated level widths (squares) found by non-linear least-squares fitting of a series of 
Lorenzians to the local density of states at the center of a triangular corral, side 220.5 A constructed 
from 72 Fe atoms on Cu(lll) (the structure reported in Ref. |22|). Calculation parameters: 
ImS/ = 0, m* = 0.38, im*t = 1, Eq = —0.44 eV. The lines are plots of the circular corral linewidth 
relationship, equation for: dashed line, radius 81.9 A; and solid line, radius 63.7 A. The inset 
shows the atomic structure of the triangular corral, and the size of the two circles. 

suits in the sparsest spectrum, which simplifies the extraction of linewidths. This point 
is discussed more fully below. Secondly, it is not generally possible to obtain an analytic 
expression for the spectral linewidths in the case of non-circular resonators. Nevertheless, 
non-circular resonators are important, for example the atomic structure of the (111) surfaces 
favours the natural formation of triangular and hexagonal resonators. 

Figure 0] shows the calculated linewidths due to lossy boundary scattering of electron 
states with amplitude at the center of a 72 Fe atom triangular corral on Cu(lll), side 
d = 220.5 A. An analytic treatment for this geometry is not possible, and so these widths 
have been obtained by fitting a series of Lorenzians to the local density of states which is 
obtained by numerically evaluating the Green function using equations I2a|2 bl The resonant 
levels occur at energies that are close to those of the hard wall limit, E p>q = 87r 2 h 2 (p 2 + q 2 + 
pq) I (3m*d 2 ) where p = 1, 2, 3, . . . , < q < p — 1, and p — q ^ 3 x integer.— Also shown in 
the figure is the linewidth relation (equation (|lip) for circular corrals with radii S = 81.9 
A, corresponding to the same area as the triangular corral, and S = 63.7 A, corresponding 
to the largest enclosed circle. In each case in evaluating (f7b|) we have used for //, the linear 
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density of corral atoms, the value from the actual triangular corral, /i = 72/ (3d). The widths 
in the same-area circular corral underestimate the triangular corral widths by some ~ 25%, 
but those for the inscribed circular corral provide a good description. Calculations on other 
hexagonal and rectangular resonators confirm this result. The lossy scattering contribution 
to the level widths at the center of non-circular resonators are approximately given by the 
appropriate circular resonator width relation (fTT|) or ()19j) . using for the radius that of the 
largest enclosed circle. Therefore when designing potential resonator structures for lifetime 
studies, lossy-boundary scattering effects can be estimated using these expressions, although 
more detailed calculations taking into account the actual geometry are needed for subsequent 
quantitative analysis of experimental spectra. 

V. THERMAL AND INSTRUMENTAL EFFECTS 

For completeness, we note that in addition to the lossy-scattering contribution, the ex- 
perimental linewidth also includes contributions due to instrumental and thermal effects. 
These must also be borne in mind when considering the use of resonators for lifetime de- 
terminations. The thermal broadening is induced by the temperature dependence of the 
Fermi-Dirac distribution function,— which enters via the requirement that electrons tunnel 
in STS from occupied states in the tip/sample to unoccupied states in the sample/tip. The 
effect on the differential conductivity is to convolute with the derivative of the Fermi-Dirac 
distribution, 



The thermal broadening function shown in Figure |3] has a FWHM of 3.5kT. This is a 
relatively small contribution at low temperatures (=1.5 mV at T = 10 K). 

The differential conductivity is normally measured using a lock-in technique, which re- 
duces phase- incoherent noise contributions.— A sinusoidal voltage modulation V m cos(ujt) is 
superimposed upon the tunneling voltage V, and the signal at frequency ui recorded: 



The measured voltage is the differential conductivity convoluted with the instrumental func- 




(20) 




(21) 
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X»(V)/X„(0) 




l.7V m 







FIG. 5: (a) Thermal broadening function, resulting from the temperature dependence of the oc- 
cupation of electron states, (b) Instrumental broadening associated with the use of a modulating 
voltage and lock-in technique for recording the differential conductivity. 



tion Xm, 



V^m - v*/*K \v\ < v m 
o \v\ > v m . 

The FWHM of this instrumental broadening function Xm> shown in Figure is 1.7V m . 



(22) 



VI. DISCUSSION 



The results of the previous sections permit the design and use of resonator structures for 
the measurement of intrinsic lifetimes. To illustrate this, we consider circular atomic corrals 
constructed from Ag adatoms on Ag(lll), for which a scattering phaseshift 5 = 0.75 + 0.42z 
has been measured.^ The basic idea is that resonators with different dimensions result in 
different series of energy levels, at energies given by (JSJ) (or more accurately by numerically 
solving (j7J)). By choosing an appropriate radius S, levels can be positioned at specific 
energies e = E n . Using STS to measure the corresponding level width T, equations (fTUj) and 
(lllj) plus knowledge of the thermal and instrumental broadening effects (section EJ may 
then be used to identify n(e) = h/Y\ n . 

Figure El compares theoretical estimates for the intrinsic many-body widths of the Ag(lll) 
Shockley surface state with the boundary-loss contributions Ten that arise in variously-sized 
circular corrals. For the smallest corral of radius S = 10 nm Tq (squares) is significantly 
greater than Tj (solid line) for most energies. In practice this will make harder the accurate 
determination of T\ and hence lifetimes 77 (= h/Ti) from STS-measured linewidths, which 
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FIG. 6: Linewidth contributions to states at the center of circular Ag adatom corrals on Ag(lll). 
Intrinsic width due to e-e and e-p scattering (solid line — compiled from theoretical values reported 
m Refs. 0QQ)- 

Lossy boundary scattering: 60 atoms corral, radius 100 A(squares); 120 atoms, 
radius 200 A(circles); 240 atoms, radius 300 A(lozenges). Inset: the LDOS normalised by pc = 
m*/irh 2 at the center of the 120 atom corral. 

will be the sum of I\ and T c , plus the effects of the convolution with xt (EHj) and Xm 
(I22j) . In the 20 nm radius corral the intrinsic and boundary-loss linewidth contributions are 
comparable for most energies, making it a better candidate structure for intrinsic lifetime 
determinations. One might consider going further, noting that for a fixed linear density 
of atoms making up the corral, fi = N/(27rS), the linewidth Tq decreases linearly with S, 
encouraging the use of corrals of even greater radius. For almost all energies the lossy- 
scattering contribution Tq in the 30 nm corral shown in Fig. El is smaller than the intrinsic 
linewidth Ft. However, from (jHJ we see that the level spacings vary as 1/S 2 , decreasing with 
radius more rapidly than Tq, so that in larger corrals the resonances merge, and the levels 
lose their integrity. This is already apparent at low energies in the 20 nm corral. The LDOS 
for this case is presented in the inset in Fig. and shows that the lowest two levels are 
only just separated. If typical thermal and instrumental broadenings are included, then the 
lowest level becomes an indiscernable shoulder. Hence lifetimes r l g are best measured in 
Ag adatom corrals with radii up to ~ 20 nm, with smaller resonators used for lifetimes at 
energies towards the bottom of the band of surface states. 

We conclude by summarising the main results of this work. We have derived analytic 
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expressions for the contribution to the spectral linewidths of electronic states measured 
using scanning tunneling spectroscopy at the centre of circular atomic corrals and adatom 
and vacancy islands that arises due to lossy boundary-scattering. Correcting measured 
linewidths for these contributions as well as thermal and instrumental broadening effects 
enables intrinsic linetimes due to electron-electron and electron-phonon scattering to be 
determined. The expressions that we have obtained are straightforward to evaluate, and 
also provide an estimate of lossy- confinement effects in non-circular resonators where the 
relevant radius is that of the largest enclosed circle. Using our linewidth expressions it is 
possible to design resonator structures appropriate for lifetime studies. 
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